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Abstract 

The time evolution of a dilute atomic boson-fermion mixture is studied at 
zero temperature, when the external magnetic field is brought close to a in- 
terspecies Feshbach resonance. On short time scales we find that before a 
collapse or phase separation of the mixture can occur, oscillations between 
fermionic atoms and heteronuclear molecules dominate the dynamics. Apply- 
ing the model presented here to a trapped mixture of 40 K- 87 Rb atoms, we 
find that the predicted oscillations can be revealed in degenerate mixtures of 
these atoms. 
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I. INTRODUCTION 



One of the most attractive features of cold atomic gases is that the interatomic scattering 
length is under very good experimental control. Using the presence of Feshbach resonances, 
the scattering length can be tuned from highly repulsive to strongly attractive by chang- 
ing the applied magnetic field. Sympathetic cooling techniques allow to cool mixtures of 
fermionic and bosonic atoms down to extremely low temperatures where the quantum me- 
chanical nature of both bosons and fermions becomes apparent, resulting in a Bose-Einstein 
condensate coexisting with a filled Fermi sea [1,2]. When the interspecies interactions are 
strongly attractive, changing the number of bosons alters the mean-field interaction strength 
and can result in a collapse of the Fermi gas as observed by Modugno et al. [3]. This insta- 
bility can also be triggered by changing the scattering length itself, making use of a Feshbach 
resonance. The divergence of the scattering length at a Feshbach resonance is due to a mag- 
netically tunable molecular state whose binding energy vanishes at resonance. This means 
that close to this magnetic field value, the molecular state can be populated. A number of 
experimental [4] and theoretical works [5-8] have been devoted to study the conversion of 
atoms between the atomic and molecular condensates and the thermal cloud. 

In this work we study the evolution of an initially noninteracting boson-fermion mixture 
at zero temperature when a sudden change (on a time scale much faster than any other 
time scale in the system) in the magnetic field brings it close to an interspecies Feshbach 
resonance. The (fermionic) molecules first form out of a condensate boson and a fermionic 
atom with an energy smaller than or equal to the Fermi energy. We will refer to such 
fermionic atoms (molecules) as 'slow' fermions (molecules); and to fermions or molecules 
with more energy than the Fermi energy as 'fast fermions or molecules, respectively. The 
decay of a molecule can then either lead back to a condensate boson with the same 'slow' 
fermion or to a thermal boson accompanied by a 'fast' fermion. Energy conservation is 
not violated because the interaction energy is changed. The latter process shows that the 
stepwise change of the magnetic field can lead to heating of the system, while the former 
demonstrates that coherent oscillations between fermionic atoms and molecules are possible, 
when a fermionic atom interacts with the condensate as a whole. Due to multiple scattering 
of fermionic atoms, bosons and molecules, the time evolution of the 'fast' part of the system 
is expected to be rather chaotic. This is why we preferred to focus on the 'slow' fermionic 
atoms and molecules. The energies involved in the formation of molecules are easily two 
orders of magnitude bigger than the Fermi energy and thus we conclude from the difference 
in phase space volume that it is rather unlikely that the 'fast' decay products of a molecule 
will come back in the Fermi sphere. This assumption allows us to study the 'slow' fermionic 
atoms and molecules subsystem. 

Our study is based on the model Hamiltonian introduced in Sec. II. We make some 
approximations that allow us to write the time dependence of 'slow' fermionic atoms and 
molecules in analytical form. The results obtained with this simplified model are discussed 
in Sec. Ill, together with the experimental conditions needed to make the predictions ob- 
servable, using data from 40 K- 87 Rb scattering. 
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II. DESCRIPTION OF THE DYNAMICS 



Experimentally, the scattering length as a function of the magnetic field close to a Fes- 
hbach resonance can be well parameterized by 



. / AB \ 

Oscatt [B) — a bg II - — - — I 



(1) 



where a bg is the background scattering length, B the position and AB the width of the 
resonance. 

A good description of the scattering physics near a Feshbach resonance can be obtained 
by explicitly introducing molecules [9,10]. These molecular states are the bound states of an 
interatomic potential in a closed channel and their energy relative to the threshold is given 
by 

u(B) = (B-B )Afji, (2) 

where A/j, — f i molcculc _ ^fermion _ ^boson j g ^ e difference in magnetic moment between the 
molecular and atomic states. Denoting the coupling strength between atoms and molecules 
by g m , the Hamiltonian for a trapped gas can be written as 



H= / d 3 x 



V} (x) H° f (x) V/ (x) + 4 (x) H° b (x) ^ b (x) + (x) H° m (x) ^ m (x) 
+9m ( x ) *l>f ( x ) A (x) + g* m ipl (x) ip\ (x) ^ m (x) , (3) 

where H° f>b (x) = -h 2 V 2 /2m f , b + V™, (x) and H° m (x) = -h 2 V 2 /2m m + V t ™ p (x) + z/ barc (B) 
are the Hamiltonians of the non-interacting gases of fermionic atoms, bosons and molecules, 
respectively. The corresponding masses of these atoms and molecules, and the trapping 
potentials confining them are indexed by /, b and m to indicate fermionic atoms, bosons 
and molecules, respectively. In the last term we stress that the parameter z/b are appearing 
in the Hamiltonian (3) is the bare molecular binding energy [10]. The true binding energy 
v (B) given in (2) differs from this quantity due to virtual dissociations of the molecules. 
Background scattering is omitted, because we assume that the system is brought close 
enough to resonance so that the resonant contribution is much more important. 

Starting from Hamiltonian (3), we wish to study the dynamics of the mixture when the 
scattering length is suddenly brought close to a Feshbach resonance. In order to grasp the 
basic physics, some approximations are made. First of all we assume that there are much 
more condensate atoms than fermions, so that the effect on the condensate can be neglected. 
Also the depletion of the condensate is presumed to be small, which allows for the Bogoliubov 
approximation. Finally we assume that the distance that a particle travels during a molecule 
formation process is small compared to the harmonic oscillator length. The latter permits us 
to restrict ourselves to homogeneous systems at first and use a local density approximation 
to take the inhomogeneity due to the confinement potential into account at a later stage. 
All these assumptions will be justified in Sec. III. 

We thus consider a homogeneous boson-fermion mixture without interactions before 
t — 0, the moment when the magnetic field is suddenly brought close to resonance and 
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the interactions between the atoms start taking place. At t ^ 0, the initial state consists 
of condensate bosons, 'slow' fermionic atoms (filling the Fermi sphere), and no molecules. 
From t — onwards, the evolution of this mixture is described by (3). The evolution will 
consist in the formation of molecules and their subsequent decay. Pauli blocking will limit 
the decay channels. However, the molecular binding energy v is, even for magnetic fields 
that are only 0.1 G away from resonance, much larger than the Fermi energy (since v = 
0.1 G XfJ'B/ks ~ 6.7// K, with hb and ks respectively the Bohr magneton and the Boltzmann 
constant, whereas E F is of the order of a few hundred nK) and blocking effects will be 
small. As a consequence, a single particle description of the fermions should be a good first 
approximation. The simplified Hamiltonian for this fermion-boson mixture reads 



H = 



I w 15^ (k) * (k) + (k) * (k) + (£ + *~ <k) *■ <k) 

77—3^ (k) i> s (k) + g m / J 3 3 ^ (k/ + k 6 ) V/ (k f ) ^ b (k b ) + h.c. (4) 
{2,71) J {Ztt) (An) 

Here, g m \p^b has the dimension of energy and we will use 

r = i x g m y/n^/H, and m = g m m r /(h 2 n b ,& ) (5) 

as our dimensionless quantities, with m r = (mj 1 + m^ 1 ) 1 the reduced mass of a fermion- 
boson pair. The first line in (4) represents the free gases of fermionic atoms, bosons 
and molecules, respectively. The second line in (4) describes the interaction converting 
a fermionic atom and a boson into a molecule ('slow' or 'fast') and vice versa. In the sin- 
gle particle picture, we should study the time evolution of each fermion in the Fermi sea 
individually. However, because the Fermi level is much smaller than the molecular binding 
energy, the initial kinetic energy of the 'slow' fermion does not influence the final result for 
the dynamics and can be neglected. Consequently all the 'slow' fermions will have a similar 
time-evolution due to the interaction with the condensate, irrespective of the specific value 
of their momentum. Thus, it will suffice to analyze the properties of a fermion with initial 
momentum k — to know the evolution of all of them. 

An appropriate wave function to describe one fermion that is initially at rest is 

\1> (*)> = [/ (0 $ (*/ = 0) |0> + m (t) {k m = 0) |0> 

+ J d%w(k f ,t)^ f (k f )4(-k f )\0)] x \N b ), (6) 

where / (t) and m (t) represent the amplitude of the bare fermionic atom and the bare 
molecule. The third term, containing the amplitude w(kf,t) arises from molecules that 
dissociate in a fast fermionic atom and a thermal boson, which is an important process. 
Higher order processes where the fast fermionic atom forms again a molecule that dissociates 
are neglected in this approach. Diagrammatically speaking, this comes down to retaining 
only the diagrams of Fig. 1 (a) and neglecting the contributions like those in Fig. 1 (b). 

The Schrodinger equation ihd t = H \ip) allows to find the time evolution of the 
amplitudes in wave function (6): 



4 



ih^f (0 = 9mV^b m (*) > ( 7 ) 

d f d 3 ~kf 

ih—m(t) = v haie + g m Vnbf (*) + 0m y ^-^w (&/,£) , (8) 

2^,— w (A;/, i) = u> (/c/, t) + g m m (t) . (9) 

We choose as initial conditions / (t = 0) = 1 and m (t — 0) = w (&/,i = 0) = 0. This cor- 
responds to the initial state wave function describing one 'slow' fermionic atom \ipinit) = 
( kf = 0) |0) . Note that y/rtb is assumed time-independent. The linear differential equa- 
tions (7)-(9) can be solved with a Laplace transform, yielding 

/ (t) = ? v e tW0T 

co \2u + v + £ (2m) 3/2 

Jo ^ 2 (co( U -i»-l) 2 + ^ 



and 



m ( T ) = \ e lWoT 

3/2 



16tt 2 



where ujq is a root of uj(uj + v + ^ (2m) ' ^/a;) — 1 = and i> is a dimensionless parameter 
related to the real molecular binding energy vby u = v/ (g m y/nb). In evaluating the integral 
in (8), we have used the relation between the real and bare molecular binding energy, v and 
z/bare, for a momentum cutoff K [10]: v = z/b ar e — Kg^rrir/ {jih) 2 . 



III. RESULTS AND DISCUSSION 

A. Atom - molecule oscillations 

It is instructive for the understanding of the time evolution to investigate the eigenstates 
of the Hamiltonian (4) describing the system in which the time evolution takes place. The 
ground state for a single fermionic atom in a system described by (4) can be investigated 
by using (6) as a variational wave function with f,m,w(kf) as variational parameters. The 
variational expression for the ground state energy is 



w."-(M) = (/ . m+ro . /)+ y^| |m(M| 



2 



gm^/Ub J (2tt) 2m 

d% 



+ J—^(rn*w(k f )+w*(k f )m), (12) 
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where kf is now measured in units n^ 3 . Minimizing this expression, with the constraint 
— 1 we fi n d f° r the ground state energy the equation 



d 3 k 



E\E-D - I — ^- w \ -1 = 0, (13) 



2m 



where E = E/(g my /n^), which has as lowest solution (relative to the bottom of the fermion 
band) E = —ujq. The ground state is a bound state representing a coherent superposition 
of the states corresponding to the fermion as an atom and as a part in a molecule: 



Abound) = A/" 



— 4>\(k f = 0)+4>l(k m = 0) (14) 



+ f — 1 —^4(k f )4(-k f ) 



10), 



with Af a normalization constant. Similar bound states are found in the photodissociation 
of a negative ion [11]. Equation (13) has only one solution for E negative. The solutions 
of (13) with E = uj > form a continuum, because we consider the limit of an infinite 
system. The eigenstates \ipexc iyj)) corresponding to these solutions are the excited states of 
the Hamiltonian (4). 

The time evolution of the initial state iV'imt) can then be written in the basis of the 
eigenstates of (4) as 

IV'(O) = Abound) e^°' (ll>bound\i>init) 

POO 

+ / \^exc(u)) e~ iuJt (ip exc (uj)\^ ini t) dcu. (15) 
Jo 

In Fig. 2, the modulus square of (ipboundl^init) , is shown for several values of v in the left 
panel. In the right panel, the modulus square of the projection of the bound state on the 
molecular state, | {Abound] V'm (^m — 0) |0) | , is shown. Note that expression (15) has the 
same structure as expressions (10), (11). This is to be expected, because the amplitudes f(t) 
and m(t) represent the projection of \ip{t)) on the 'slow' fermionic atom component and the 
'slow' molecule component: 

f(t) = (0| i> s (k f = 0) |^)>, 
m(t) = (0\^ m (k f = 0)mt))- 

The two terms appearing in \ip(t)) , expression (15) give rise to the corresponding two terms 
in f(t), expression (10), and in m(t), expression (11). The first term (the first line in 
expressions (10), (11)) is related to the bound state and the second term (containing the 
integration over uj) is related to the excited states. Because of the continuous character of 
the excited states as a function of u>, the terms related to the excited states vanish for large 
times. 

We find that for large the conversion from 'slow' to 'fast' fermions is greatly reduced. 
This can be understood from the fact that for large \u\, (ipboundl^init) becomes very small 
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(as can be seen in Fig. 2) and the initial state coincides nearly with an eigenstate of the 
system, \ipi n u) ~ \ip exc (ui = 0)) so that \ip{t)) ~ \ipmit) and fit) remains close to 1 (its initial 
value) . 

The fraction of 'slow' molecules or 'slow' fermionic atoms, obtained from (10) and (11) 
respectively, displays damped oscillations resulting from the interference between the terms 
related to the bound state and the terms related to the excited states. Fig. 3 shows the 
evolution of the densities of 'slow' fermionic atoms and 'slow' molecules for different values 
of the molecular binding energy u and mass fh. For the chosen values of the mass, the 
oscillations are visible for negative molecular binding energy v < 0, while for u positive they 
are strongly damped. 

We observed that the larger rh, the fewer molecules are formed. Indeed, according to 
Fermi's golden rule, a bigger reduced mass means a faster dissociation rate (~ vr?J 2 ) of the of 
molecules. Remarkably, also for m ^> 1, the fraction of 'slow' fermions reduces substantially 
over time even though at all times few molecules that can dissociate to 'fast' fermionic 
atoms are present. This effect can be understood since, even though few molecules are 
present at any time, they form and dissociate quickly. This raises the question whether our 
approximate wave function (6) is good enough. For low m, it is obvious that the dissociation 
of molecules will be slow and strong Rabi oscillations between fermionic atoms and molecules 
occur. The third term of (6), containing the factor w(kf,t), is then relatively unimportant 
and corrections on it will not affect the result a lot. For high m, the molecules dissociate 
rapidly and the contribution to \ip (t)) coming from the pairs of 'fast' fermionic atoms and 
bosons (last term from equation (6)) must be big. But in this case the corrections to the 
contribution of these pairs -through new formation of fast molecules out of a fast fermionic 
atom- are small, because few of those molecules are formed. Diagrammatically, this means 
that diagrams such as those in Fig. 1 (b) are negligible. Then for intermediate rh one can 
also expect (6) to be a good approximation, at least in order to calculate the fraction of 
'slow' fermions. 



B. Relation to experiments on °K- Rb mixtures 

To check whether the picture described here, with fermions that form a bound state with 
the condensate, is correct, one has to be able to measure preferably the fraction of 'slow' 
fermionic atoms. This should, after some transient behavior, go to a constant value and not 
to zero. They can be detected by a time of flight expansion after a magnetic pulse has been 
applied. 

In 40 K- 87 Rb plenty of Feshbach resonances have been predicted [12,13] allowing us to 
estimate quantitatively the magnitude of the model parameters. Therefore, we make use 
of g m = ^/2irh 2 ab g AB Afi/m r , with a bg the background scattering length. As an ex- 
ample we take the resonance between the (9/2,-9/2) + (1,1) hyperfine states around 
741.2 G- For n b = 10 14 cm~ 3 , the frequency g m ^fnbjh is equal to 7.72 x 10 5 Hz, which is 
much larger than the harmonic oscillator frequency. Thus we can expect that the con- 
finement is relatively unimportant for the atom-molecule oscillations. We find in this case 
v = 33.13 x {B - B ) (G) /y/n b (10 14 cm~ 3 ) and further m = 15.47. The values used for the 
parameters of Fig. 3 are for this resonance with n b = 10 14 cm -3 and B — B = ±0.1 G on the 
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left, and on the right for the (9/2, —7/2) + (1, 1) hyperfine states with a Feshbach resonance 
around 522.9 G- 

The current results have been derived in the framework of a translation invariant system 
whereas in the experimental setup the inhomogeneity of the trapped gasses is apparent. The 
time evolution has been calculated for homogeneous mixtures, and it is necessary to estimate 
what the effects due to the confinement are. Considering the case with much more bosons 
than fermions, we can neglect in a first approximation the change of the boson density. Then 
one can find in a mean-field approximation the redistribution speed of the fermions in the 
trap. The local density approximation is applicable only if during the time scale of a Rabi 
oscillation between the fermionic atom state and the molecular state, the fermions have not 
moved over a substantial fraction of the trap length. In the Thomas-Fermi approximation for 
the bosons, the mean field potential felt by a fermion is 27rh 2 a scatt n^/m r x (1 — r 2 /R 2 -,) , with 
n° the condensate density in the origin and Rq the radius of the condensate. This results 
in an additional confinement potential (inverted for a scatt > 0) acting on the fermions. For 
the resonance used above, for n° = 10 14 cm~ 3 and R c = 10 /im and making use of (1), its 
frequency equals uobf = 4.95 x 10 3 s _1 /-\/£. For v not too small, this is much lower than 
the frequency scale of the molecule formation j m y^//i = 7.72 x 10 5 Hz, which means that 
the dynamics due to molecule formation is fast compared to the mean-field collapse, and 
the local density approximation can be used. For v — 1, an atom migrates 10% of the trap 
length in 0.75 ms whereas the period of a Rabi oscillation is around 1.3 /is. 

Within the approximations described above, we plot in Fig. 4 the density of slow fermions 
as a function of time and distance from the center of the trapping potential in the case that 
the fermion cloud has a radius 1.5 times bigger than the condensate. To use the formalism 
presented above, space was subdivided in discrete boxes in which the gases were presumed 
homogeneous, as in the local density approximation. As discussed above, no redistribution 
of fermions across volume elements is considered so that each volume element can be treated 
independently of the others. In the lower panels, the total fraction, averaged over space, of 
'slow' fermionic atoms is compared to the local fraction of 'slow' fermionic atoms in the center 
of the trap. One can see that the inhomogeneity decreases the amplitude of the oscillations 
in the total fraction of 'slow' fermionic atoms. The reason is that the condensate is the 
most dense in the middle, which increases the efficiency of the atom-molecule conversion, 
and also because r, fh and v depend on rib ( r ) , so that on different places in the condensate 
oscillations with different frequencies occur. 

IV. CONCLUSIONS 

We have made a study of the dynamics of a fermion-boson mixture resulting from a 
fast switch in the magnetic field, bringing it close to an interspecies Feshbach resonance. 
Our analysis shows that damped oscillations between fermionic atoms and molecules will 
occur on the typical time scale of h/(g m ^fnb) (~ 1/xs for a 40 K- 87 Rb mixture). Such a type 
of oscillations have recently been observed in single component BEC's near a Feshbach 
resonance [4]. The decay of the molecules in pairs of 'fast' fermionic atoms and bosons leads 
to heating of the system, but a core of 'slow' fermions will remain, due to the existence 
of a bound state of the fermions coupled to the condensate and the continuum of thermal 
bosons. Using the latest data on resonances in 40 K- 87 Rb scattering [12,13], we find that 
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the predicted oscillations between slow fermions and slow molecules should be observable 
both in a homogeneous and inhomogeneous system. The redistribution of the densities are 
expected to occur much slower than these oscillations. 
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FIGURE CAPTIONS 

FIG. 1: Feynman diagrams that describe the correction to the Green function of the 
heteronuclear molecules in the mixture. The terms in the left part (a) are taken into account 
in the theory of this paper. Multiple dissociations as depicted in the right part (b) are 
neglected. 

FIG. 2: Number of molecules (right) and fermionic atoms (left) for t — > oo as a function 
of the molecular binding energy v for different values of the mass m. 

FIG. 3: Time evolution of the densities of 'slow' fermionic atoms and molecules in 
a homogeneous system for two different values of the molecular binding energy and the 
mass. The molecular binding energies correspond to detunings of the magnetic field of 
B—Bq = ±0.1 G- For rif, = 10 14 cm -3 the dimensionless detunings and the masses correspond 
then to the predicted 40 K- 87 Rb resonances in the (9/2,-9/2) + (1,1) hyperfinestates at 
741.2 G (left) and to the (9/2, -7/2) + (1, 1) at 522.9 G (right). 

FIG. 4: Time evolution of the density of 'slow' fermionic atoms when a noninteracting 
boson-fermion mixture in a harmonic trap is suddenly brought close to a Feshbach resonance. 
In the upper panels, the density of 'slow' fermionic atoms is shown as a function of the 
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distance to the center of the trap and in the lower panels, the spatially averaged fraction 
of 'slow' fermionic atoms is compared with this fraction in the origin. The inhomogeneity 
decreases the spatially averaged oscillations because the amplitude and frequency depend 
on the condensate density and consequently vary in space. The parameters used are the 
same as for the right panels in Fig. 3, m — 3.8 and v = ±2.5, and correspond to a realistic 
situation in 40 K- 87 Rb mixtures. 
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